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Abstract 

The Darboux-Dressing Transformations are applied to the Lax pair associated to the 
system of nonlinear equations describing the resonant interaction of three waves in 1+1 di- 
mensions. We display explicit solutions featuring localized waves whose profile vanishes at the 
spacial boundary \x\ — oo, and which are not pure soliton solutions. These solutions depend 
on an arbitrary function and allow to deal with collisions of waves with various profiles. 
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1 Introduction 



The propagation in 1+1 dimensions of three resonating waves is modelled by the following system 
of three coupled PDEs 

DnXn =5nXn+lXn+2> fl = 1, 2, 3 mod (3) , (1) 

where the three dependent variables x„ = Xn(x,t) are generally complex functions of two real 
independent variables, x, (space) and t (time). The three differential operators D n , n = 1, 2, 3, in 
their linear (i.e. left-hand-side) part are first-order and read 

D n ■= d t + c n d x , (2) 

where the three real constants c n are taken all different among themselves. Each c n is the char- 
acteristic velocity at which the field Xn propagates in the space-time region where no interaction 
takes place. Here and hereafter a subscripted independent variable denotes partial differentiation 
with respect to it, asterisks denote complex conjugation and all indices are considered modulo 
3. The three constants g n are generally complex and have the significance of coupling constants. 
Their values depend of course on the physical context to which the 3Wave Resonant Interaction 
(3WRI) equation applies. It is remarkable however that, if these coupling constants satisfy certain 
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conditions (see below), then the 3WRI equation is integrable. In this paper we consider only 
such case and apply Darboux (or "dressing" ) transformations to construct explicit solutions of the 
system (Q. 

The conditions on the complex constants g n which are necessary and sufficient for the integrability 
of the 3WRI equation Q are: 

g n = <r n \g n \ exp(ir) , cr£ = 1 , n= 1,2,3. (3) 

Here the three signs a n and the three moduli \g n \ are arbitrary and the condition is that the (also 
arbitrary) phase T is common to the three g n 's. It is worth pointing out here that these conditions 
do single out only one integrable system in the following sense: the structure of the system JJJ is 
left unchanged by the transformation 

Xn(x,t) = a n Xn(x,t) , x = Ax + Bt . i — Cx + Dt , (4) 

where the three arbitrary constants a„ are complex and the four constants A, B, C, D are real and 
arbitrary except for the requirement AD — BC ^ 0; indeed it is clear that such transformation 
only changes the value of the six constants c„ and g n . It is easily seen that the transformation (@J, 
as it should, takes g n 's satisfying the integrability conditions © into different g n 's which satisfy 
the same conditions, and that any set of values of g n 's satisfying © can be transformed into any 
other set which satisfy ©, for instance g\ = g<z = g^ = 1. Without any loss of generality, hereafter 
we choose the coupling constants to have the expressions 

9n = Cn+i - c n +2 , n = 1 , 2, 3 mod (3) (5) 

in terms of the velocities c„ with the implication that gi + g% + gs = and, therefore, that the g n 's 
cannot have all the same sign. The reason for this choice is practical as the Lax pair associated to 
the system pQ in this case takes a simple look and it reads 

$ x = (-ikB + U(x,t))V, (6a) 



y t = (ikA + V(x,t))V, (6b) 
where A and B are two diagonal, real, traceless matrices 

A = diag{ai, a 2 , a 3 } , B = diag{6i, b 2 , b 3 } , 

(7) 

A = A*, B = B* , tr(A) = tr(B) = Q . 

The 3x3 matrix ^> = ^(x,t,k) is the common solution of the two ODEs © while U(x,t) and 
V(x, t) are off-diagonal matrices whose entries are related to the three wave fields Xn{ x i t)'- 

( X3 -X*2 \ ( -C3X3 C 2 X*2 \ 

U= \ -Xt xi , V = c 3X l -cixi • (8) 
\ X2 ~Xi / \ -C2X2 C1X1 / 

Moreover the three characteristic velocities c„ which appear in are related to the six real 
constants a n and b n by the relations 

- - a »+ 1 " a "+ 2 n = l,2,3 mod (3). (9) 



b n +i - b 



n+2 



The integrable character of the 3WRI equation QJ, together with JSJ, follows from the compatibility 
(^ xt = ^>tx) of the two linear equations ©■ 

The system as it models the interaction of three resonating waves, has attracted much attention 
in the study of nonlinear waves (see for instance [2])- It can be derived in a quite general way 
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via multiscale analysis of a large class of nonlinear wave equations with weak dispersion and 
nonlinearity (see for instance [3]). As such, it represents a universal model which is able to capture 
the lowest order corrections, due to nonlinear effects, to linear propagation. It is here appropriate 
to point out two special properties of the equation which make this model peculiar with respect 
to most of the intcgrablc nonlinear wave equations (f.i. Nonlinear Schrodinger equation, Korteweg- 
de Vries equation, Sine-Gordon equation and others). The first one is that the system JJJ has no 
(linear) dispersion (see J2J). The second property is that no self-interaction occurs, namely, the 
forcing of each wave (i.e. the nonhomogeneous term in the propagation equation |JT)|) is just the 
product of the (complex conjugated) amplitudes of the other two waves, see the right hand side 
of JTJ. If we assume that the spatial profile of each one of the three waves is well localized in the 
far past (i.e. as t — > — oo), then these two properties imply that each wave profile Xn translates 
with no deformation with its characteristic velocity c„ up to the time it encounters an other wave. 
Then, for some time, the interaction takes place where two waves (or the three of them) have 
spatial overlapping, but eventually, in the far future (i.e. as t — ► +00), each wave is again free to 
move with no deformation with its own characteristic velocity. Generally, the three wave profiles 
at t = —00 arc different from those at t = +00 since the interaction at intermediate time causes a 
change of the profile. In the (x, t) plane the asymptotic behaviour is therefore expressed by 

Hrn Xn(y + ct,t) = 0, c^c n , lim Xn{y + c n t, t) = xl ±} (y) , (10) 

t — >±oo t — *±oo 

where y is the spatial coordinate in the moving reference frame, and the functions Xn ~ c nt) 
are the asymptotic field profiles. Of course, the relevant physical problem is computing the three 
wave packets in the future, Xn(y), once the packets Xn(y) m the past are given. This large 
time asymptotic behaviour of the three waves may be referred to as "asymptotic freedom". In 
terms of the spectral data associated to the three functions Xn{x,t) via the spectral problem 
^ x = (—ikB + U(x. t))^, see (|6all . the spectral components of the asymptotic states Xn^iv) ma y 
well be on the continuum spectrum (wave-packets) other than on the discrete spectrum (solitons) . 
These phenomena are quite in contrast with the well-known "separation" process occurring for 
dispersive nonlinear waves where only the discrete spectrum components (solitons) show up as 
t —> ±00, since the continuum spectrum components (wave packets) separate from solitons, disperse 
away and vanish (see for instance [4]). The asymptotic freedom behaviour of the three interacting 
waves may be given a spectral formulation [2], which is however not discussed here. 

In Section 2 we collect formulae which provide us with the tool of the Darboux-Dressing Transfor- 
mation (DDT), while in Section 3 we obtain, via the DDT, explicit solutions of the 3WRI equations 
Q). In order to show the asymptotic freedom for wave-packets, we deal with explicit solutions 
which are not pure soliton solutions, this being the main result of this paper. An outlook with few 
remarks are the content of Section 4. 



2 The Darboux-Dressing Transformation 

Let us now turn our attention to the method of construction of explicit solutions of the 3WRI 
equation Q). This construction is based on the Darboux transformation of ODEs and since this 
transformation is well-known [5] we limit our presentation below to a collection of formulae and 
propositions. To this purpose it is convenient to observe first that, without any loss of generality, 
we can give the constants b n and a n , n = 1,2, 3 which appear in the Lax pair of equations (jfiaj) 
and, respectively, l|6bj) . see also JjJ, the following expressions 

bn = g (ffn+i - 9n+2) , a n = -{c n+1 g n+ i - c n+2 g n +2) , n = 1, 2, 3 , mod (3) . (11) 

The validity of this statement follows from the expressions © of the characteristic velocities and 
(JSJ) of the coupling constants, from the relation b\ g\ + 62 <?2 + ^3 <?3 = 0, which is implied by 
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them, and from rescaling the spectral parameter k. As a result, the spectral parameter k has now 
the dimension of wave-number /velocity. It should be also noticed, for future reference, that the 
relations can be inverted to yield the expressions 

9n = -(b n +i - b n+2 ) , n= 1,2,3, mod (3) , (12) 

which can be used as alternative to (JjjJ with the implication that b n — —c n + (ci + ci + c^)/2>. We 
then note that both matrices U and V in the Lax pair JBJ satisfy the reduction condition 

U = -U\ V = -V\ (13) 

where the dagger denotes hermitian conjugation. This property Ijl3(l induces on the solution ^> of 
the two equations © the corresponding condition 

¥(x,t,k*)y(x,t,k) = C(k,k*), (14) 

where the matrix C(k,k*) is constant, namely x- and i-independent; its value depends of course 
only on the arbitrary value V&(a;o, to, k) that the solution '5 takes at a given point (xo,to) of the 
{x, t) plane. 

Consider now a second pair of off-diagonal matrices (x, t) and (x, t) which satisfy the same 
skew-Hermitian condition (C/^) t = -C/ (0) , (y (0) )t = -V (0) as U and V (see (O); their entries, 
with obvious notation (see (JHJ), identify the three complex functions Xn\x,t), n = 1,2,3. Let 
'J'- -' (x, t, k) be a corresponding nonsingular (i.e. with nonvanishing determinant) matrix solution 
of the Lax pair ©, i.e. 

*(?) = (-ikB + [/(°))*(°) , *f } = (ikA + y(°>)tf(°) (15) 

which have the same form © with U and V replaced, respectively, by and . Assume that 
the initial condition ty(°\xo,to, k) is so chosen that the constant matrix C^°\k, k*) (see Q140 ). 

C (0) (k,k*) = ¥ 0) \x,t,k*)¥ 0) (x,t,k) (16) 

takes the same value of C(k,k*), i.e. C^°\k,k*) = C(k,k*). If both compatibility conditions, 
^ X °J = ^f[^ and ^ x t = ^tx, are satisfied, then ^ and Xn, n — L 2, 3, are two different solutions 
of the same 3WRI equation JIJ , and the matrix 

D(x,t,k) = ^{x,t 1 k)(¥ 0) (x 1 t,k))- 1 (17) 

satisfies the pair of differential equations 

D x =ik[D,B] + UD-DU {0) , (18a) 

D t = -ik[D,A] +VD-DV {0) . (18b) 

Moreover, as a consequence of the reduction conditions (|14|l and (|16|l with C = , the matrix 
D(j;, t, k) satisfies also the algebraic equation 

D^(x, t, k*) D(x, t, k) = I. (19) 

The proof of these statements is straight. 

The definition (|17|l can be viewed as a transformation of ^^°> into 

V(x, t, k) = D{x, t, fc)* (0) (x, t, k) , (20) 

which consequently yields a transformation of Xn^ (the bare solution) into Xn (the dressed solution). 
Clearly the dressing approach requires in the first place that Xn^ (x, t) and {x, t, k) be explicitly 
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known. The next step is the construction of the transformation matrix D(x, t, k) via the integration 
of the ODEs (JTSJ. This task however is not straight since the coefficients of these differential 
equations depend on the unknown matrices U and V, see (II8a|l and (|18bl) . The way to overcome 
this difficulty goes through the a priori assignement of the dependence of the transformation matrix 
D(x, t, k) on the spectral variable k. 

The simplest instance of this strategy is illustrated by the following proposition (whose proof is a 
mere exercise): if D is ^-independent, D% — 0, then D is diagonal, [A, D] = 0, [B,D] = 0, and is 
i-and i-independent, i.e. D x = D t = 0. The associated transformation of C/(°' is then a gauge 
transformation, namely D = G, U = GU^G^ 1 . 

In the following we will consider the larger set of fc-dependent matrices D(k) which (i) have a 
rational dependence on the complex variable k and (ii) have nonvanishing limit as k — > oo . As 
we consider here only rational dependence on k which can be factorized as product of simple-pole 
terms, we need to deal only with matrices D(x, t, k) which take the following one-pole expression 

D(x,t,k) = l + ^-, (21) 
k — a 

the matrix R(x,t) being the residue at the pole k — a. It is plain that the transformation 
matrices D(k) whose k — > oo limit is a nonsingular, but not necessarily unit, matrix may be 
obtained by multiplying the expression (|2I(I by a gauge transformation G. The transformation 
characterized by the matrix H21(l has received considerable attention in the literature 0; we refer 
to it as to Darboux-Dressing Transformation (DDT). Its existence in our setting is proved below 
by construction. 

In general, the way to obtain an explicit expression of the residue matrix R(x, t) depends on 
whether the pole a is on the real axis, a = a* , or not. However, in the case under investigation, 
it is easily seen that the condition a = a* implies R(x,t) — 0. Thus we need to consider only 
the case in which a is not real, a ^= a*. The starting point is the requirement that the matrix 
D(x,t,k) satisfies the algebraic condition l|I9() and the differential equations (|f 8f) . The algebraic 
condition (|f 91) entails the equation 

^+^^ = 0, (22) 

a — a* 

whose solution is 

R(x,t) = {a - a*)P(x,t) , (23) 
where the matrix P(x, t) is an Hermitian projector 

P 2 = P , P f = P . (24) 

As for the two differential equations l|I8|) . replacing D(x,t,k) with its expression 

D(x,t,k) = l+(j^p\P(x,t), (25) 

(see (12 H and (J2SJ) yields two algebraic and two differential equations. The algebraic relations read 

U = U {0) +i(a-a*)[B,P], (26) 

V = V W - i(a - a*)[A , P] , (27) 

and give the dressed quantities Xn(x,t), in terms of the bare ones Xn\x,t) n = 1,2,3, and the 
projector P (see below). The two differential equations are 

P x = (1 — P){-ia*B + U {0) )P - P(-iaB + U (Q) ){1 - P) , (28a) 
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Pt = (1 - P)(ia*A + V (0) )P - P{iaA + V (0) )(l - P) . (28b) 

Let v be an eigenvector of P with unit eigenvalue, and differentiate with respect to x the eigenvalue 
equation 

p v = v , v = I v 2 I . (29) 




By replacing then P x with the right-hand side of Q28a[l . one obtains the equation 

(1 -P)[v x - {-ia*B + U (0) )v] = 0, (30) 

which implies that the vector v x — (—ia*B + U^°')v belongs to the subspace on which P projects. 
At this point we note that, since the dimension of P is 3, we may well assume, without any loss 
of generality, that this subspace be one-dimensional. Indeed, it is easy to show that if P projects 
on a two-dimensional subspace, then I P projects on a one-dimensional subspace and the only 
change in the final result is the value of the parameter a, which goes into a*. Therefore, P projects 
on the one-dimensional subspace of the vector v: 

{X,t) ~ vi(x,t)v(x,t) ' 16 ' 

where in this notation we treat vectors as one-column rectangular matrices. This implies, see (|30|) . 
that the vector v x — (—ia*B + U^°')v is proportional to v, i.e. 

v x -(-ia*B + U^)v = hv. (32) 

On the other hand, since the vector v is identified only modulo a scalar factor, which may well 
be a function of x and t, we may choose this factor in such a way that v satisfies this differential 
equation l|32() with h = 0. 



The differential equation l|28b|) with respect to the variable t can be handled in a similar way. One 
obtains therefore that the vector v(x,t) satisfies the differential equations 

v x = {-ia*B + U [a) )v , (33a) 



A + V (0) )v. (33b) 

Once the two equations (|33() are solved, the DDT transformation matrix D(x,t, k) is finally given 
explicitly by with IpHJl. 

At this point, we conclude that the method of construction of a novel solution Xn{x, t) of the 3WRI 
equation starting from the knowledge of a given solution Xn \ x 7~t) is explicitly given by (|26|l 
with H31f) where the vector v(x,t) is given by 

v(x,t) = ^ {0) (x,t,a*)v Q . (34) 

Here (x, t, a*) is assumed to be a known solution ty( '(x, t, k) of the Lax pair (|15|) . for k = a* , 
while v is an arbitrary constant vector, 




«o = 72 , (35) 



where 71, 72, 73 are three (arbitrary) complex parameters. With these specifications (see also (|12(l 'l. 
the DDT transformation (|26|l takes the more explicit form 

Xn = Xn ] ~ i( a ~ "*)5ni — , 2 ,"| ^ i — Tn , n = 1, 2, 3 , mod (3) . (36) 
\vi r + \V2 r + mr 
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3 Solutions 



In this section we apply the DDT to known solutions x n (%■> t) of the 3WRI equation (JTJ to 
obtain other solutions of the same equation. Let us consider first the DDT within the class of 
solutions Xn(x,t) of the 3WRI system which are well localized at all times. This class is 
left invariant by the DDT if the functions Xn(x,t) vanish sufficiently fast for large \x\ at fixed 
t, linxr^ioo Xn{x, t) = (for instance, if they are in L\(R) as functions of the x variable). The 
simplest choice of the bare solution Xn^ {x, t) in this class is, of course, the vanishing one, i.e. 
Xn^ (x, t) = 0. The corresponding solution which is obtained by applying to it the DDT described 
in the previous section is the well-known one-soliton solution (see below). The more general class 
of solutions we compute here originates from the simple observation that even asking that only 
two of the three fields XI1X2 and X3 are vanishing yields a solution, which depends moreover on 
an arbitrary complex function of one real variable. Since the non vanishing field can be any one of 
the three, we introduce three different bare solutions as distinguished by the index j: 

X$ d \x,t) = 8 jn f(x n ), J = 1,2,3, (37) 

where f(y) is an arbitrary complex function of the real variable y and the three variables x n , 
n = 1,2,3, are the characteristic coordinates 

x n = x — c n t, n = 1, 2, 3 . (38) 

Remark. The distinction of the three solutions l|37|) by the index j makes sense if the labeling of 
the three fields Xn by the index n is fixed. From now on, without any loss of generality, we set 
such labeling by ordering the three characteristic velocities in the following way: 

ci < c 2 < c 3 . (39) 

Indeed we remind the reader that all other relevant parameters, namely the coupling constants g n , 
see ©, and the constants a n and b n , see which appear in the Lax equations, see and 

(|6b|) with JTJl, are functions of the c n 's. 



The first step in applying the DDT to the seed solution 1)37(1 is of course solving the Lax pair of 
equations 1(1 5|1 for each one of the three cases (with obvious notation, and see©) 



) 
= ( f( Xl ) 

) -f*( Xl ) 

f(x 3 ) 
£,(0,3) = I _ r(x3) o o 





-f*(x 2 ) 
= I 

f(X2) 



V (o,j) = _ c .f/(o,j) } 



(40) 



which correspond to the three bare solutions JHTJ. Because of the particular structure of the 
matrices U^- '^, see (I40|l . and of their dependence on x and t through the characteristic coordinates, 
see (|38|l ■ the matrix solution ^(°^) (x, t, k) takes the form 



^ J \x,t,k) = ^\xj,k) exp 



k 

i — (bjx — Oji)(l — iPj) 



(41) 



where the three matrices Pj j = 1 , 2 , 3, are the three diagonal projectors 

= I (I I . P; = I I 1 • P:; 






(42) 
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Solving the equations of the Lax pair now reduces to solving only the 3x3 matrix ODE 



- gjk (P j+1 -P j+2 ) + U^\y) 



(43) 



which is however easily seen to be equivalent to the standard Zakharov-Shabat system of two 
equations 

<hv = \Mi + f(y)h , fay = -\\<h - r{y)<j> x (44) 

for the two functions 0i(y, A) and (j> 2 {y, A). Indeed, once a solution of the system (|44|) is known, 
the 3x3 matrices <S>(i\y, k), for j = 1 , 2 , 3 , which solve (|43|l are also known and read 
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$«(y,fe)= ( o Mv.gik) -4>* 2 (y,gik*) 

cj) 2 {y,gik) 4>\(y,gik*) 



(45a) 



, 4>2{y,g 2 k) -(/}l(y,g 2 k*) 

$ (2) (y,fc) = ( 1 

4>i(y,g2k) (f> 2 (y,g 2 k*) 

. 0i(y,53fc) ~<t> 2 (y,93k*) o 

$W( y ,k) = \ My, 93k) ri(y,93k*) 
1 



(45b) 



(45c) 



For future reference, we display also the asymptotic behaviour of the solution of the Zakharov- 
Shabat system ijHj) 



0i (y, A) — -* zJ ±) (A)exp ( l -\y 



^2(2/, A) z^ ±] (A)exp ( --Ay 



(46) 



which introduces the four functions z x {\) and z^^A), two of which can of course be arbitrarily 
given while the other two can be obtained by integrating the ODEs (|44() . 



Let us now consider the dressed solution by applying the DDT given in the previous section to the 
bare solution i|T7|) . Its expression (see (HSU) is 



v%{x,t)v^ 2 {x,t) 



\v[ 3 \x,t)\2 + \vif\x,t)\2 + \vi3\x,tW 



(47) 



where the parameter 77 is the imaginary part of the complex pole a which has been introduced in 
the DDT formula (J2SJ), namely 

a — p + ir] . (48) 

According to the definition (|34|) . the three 3-vectors v^(x,t), with components u„ (x,t), are 
provided by the formula fseel|41|l) 



exp 



-a>jt)(l -3Pj) 



(49) 



where the arbitrary complex constant vector vq introduces in the solution Xn {%■> t) the three 
arbitrary parameters 71 , 72 and 73, see l|33|) . This new solution (|4TJl of the 3WRI equation depends 
on the function f(y) which characterizes the bare solution 13 7|) and its behavior may be roughly 
described as the interaction of the bare j'-th wave f(x — Cjt) with soliton-type bumps in the other 
two waves whose speed may be larger or smaller than Cj , a behavior which may be quite complicate 
at finite intermediate times. Here we do not dwell in detailing such behavior at finite times, rather 
we compute the asymptotic states as t — > ±00 whose profile and properties really matter in an 
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applicative context. According to the general definition 11U|) . we now compute, for each value of 
j, the six functions Xn^^iy) by performing the limits 

lirn xiP(y + Cnt,t)= X ^ ± \y) . (50) 

t—>±oo 

This computation is elementary but lengthy, and it is omitted. The relevant results take the 
following expressions: 



W (.. _l „ .+ a - tu.\ j_ 9„„ . "H-'itri^M 

k^xWP+h^fe)! 2 ' (51a) 



Hmt^oo Xj ; (y + Cjt,t) = f{y) + ^V9jjjj) 



\imt^-s 3 oox ( j } (y + C]t,t) = f(y) 



lim t ^ Sj+l00 Xj+i(y + Cj+it,t) 



2r]g j+1 exp{ipg j+1 y) 



Jj U 3 + 2 



lirn* 



.t->-« j+1 oo Xj + i(y + c j+ it,t) = , 



exp(-77g i+ u/) + |u i+2 



! cxp(?yg J + il/) 



(51b) 



lim t ^ Sj+200 Xj+ 2 (y + c j+2*, t) 



2 Wj+2 exp(i raj+2 j/)— — (J)( 

l7iP ex P('7Sj+2y)+|« i+1 r cxp(-7)gj +2 y) 



(51c) 



Yim t -,-s j+2 oo Xj+ 2 (y + c 3+2t, t) = . 
The symbols we have introduced here are: the three signs 

sj = sign(?75 i ) , s| = 1 , j = 1 , 2 , 3 , 



(52) 



the real parameter p which is the real part of the complex pole a, see (|48|l . the three y— dependent 
vectors (see (|45al I45bl I45c|) and 



w ij \y) = a*)i> 



and the three constant vectors 



,(i)(±) 



7i 



,(2)(±) 



,(3)(±) 



/ 

72^i a ) - j 3 z 2 (51 a) 
\ 724 ±) (5i"*) + "/ 3 z ( { h) *(g 1 a) 

' llZ^ (92a*) ~ l3z{ ±] * (92a) 
72 

\ liz[ ± \g2a*) +"f 3 z ( 2 ±) *(g 2 a) 

/ (±)/ *\ (±)*/ \ 

/ 7i^i (53a ) - 724 (53a) 

liz 2 '{g 3 a ) +J2Z{ (g 3 a) 

V 73 



(53) 



(54a) 



(54b) 



(54c) 



where the four parameters z± (gjOt), %2 (9j a ) an d z i {9j a *)i z 2 (9j a *)i f° r 3 — 1j 2, 3, are the 
values that the two functions z[ (A) and (A), introduced via the asymptotic behavior l|46l) . 
take at A = 5^a and A = gja*, respectively. 

At this point we are in the position to discuss the asymptotic in- and out-states corresponding 
to our solution. Since we aim to attract the attention of the reader on potential applications, we 
chose to have the bare wave f(x — Cjt) in the initial state (say as t — > —00). Therefore, as indicated 
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by the limit formula i|51al) . we have to set Sj = 1 for each j, and this implies, according to the 
definition l|52l) and to the signs of the coupling constants, see © and (|39[) . 



signal 



-1 , sign(5 2 ) = 1 , sign(g 3 ) = -1 



(55) 



that we pick the pole a of the DDT in the lower half of the complex plane, r] < 0, for j = 1 and 
j = 3 and in the upper half of the complex plane, r\ > 0, for j = 2. By taking into account that 
these choices of the sign of the parameter r\ entail that 



s 2 = -1 , s 3 = 1 ■ for j = 1 , 
si = -1 , s 3 = -1 , for j = 2 
si = 1 , s 2 = -1 , for j = 3 , 



(56) 



the final expressions of the asymptotic profiles of the three waves are readily read out of the 
formulae (|51a(l . I)51bjl and (|51c|l . and are 



,(2)(- 



/ (3)(- 



\x,t) 



f(x - at) 

S ( 2 1) (x-c 2 t) 


/ S[ 2 \x- Cl t) 

f(x - c 2 t) 
\ Si 2) (x-c 3 t) 



Si 3) (x-c 2 t) 
f(x - c 3 t) 



X 



(i)(+) 



X 



(2)(+) 



X 



(3)(+) 



(x,t) 



(x,t) 



(x,t) 



F^(x-dt) 


4 l \x-cat) 


F i - 2 \x-c 2 t) 


S[ 3 \x-dt) 


F( 3 ) (a; - c 3 i) 



(57a) 



(57b) 



(57c) 



In these formulae we have expressed the asymptotic profiles by means of the following expressions, 
which again are just implied by QSlafl . (|51b|) and i|51c[l : 



F^(y) = f(y) + 2 m 

S 2 \y) = 2f 1 g 2 exp(ipg 2 y) 
S 3 (y) = 2f 1 g 3 exp(ipg 3 y) 

(2) 

5*1 (v) = 27751 exp(ipg!y) 

(2) 

s 3 (v) = 2 V93exp(ipg 3 y) 

51 (y) = 2 V9i exp(jpffiy) 

5 2 (y) = 2V92exp(ipg 2 y) 



^UyW + ^Uy) 



Til 2 


cxp(-r(g2l/) + |"3 1K " 

ti4 1)(+> * 


I 2 exp(?7g2l/) 






Til 2 


exp(»)3 3 y) + |«^ 1)<+) 

72 4 2)(+) * 


2 exp(— ?7flr 3 j/) 


T2| 2 


cxp(7(giy) + |?4 


2 exp(— ?73iy) 


T2| 2 


cxp(-r)g 3 y) + \u[ 2>{ 


' | 2 cxp(r/g 3 y) 


T3| 2 


c^p(-V9iy) + \un 3K ~ 


' I 2 exp(j7gij/) 


T3| 2 


cxp(j)g 2 y) + |«i 3)<+) 


2 cxp(-77g 2 y) 



(58) 



(59a) 



(59b) 



(59c) 



Their meaning is quite transparent: the function F^\y), which always appears in an out-state, 
shows how the profile of the incoming wave-packet f{y) changes because of the collision with 
solitonic waves, while the functions Sffl (y) show the standard sech-shape of the solitons (see below) 
appearing in the incoming and/or outgoing states. The three different collision processes described 
by these asymptotic states, (|57a|) . (|57bfl . (|57c|) . are controlled by various parameters introduced 
via the DDT, namely the complex pole a = p + ir/ and the complex components jj of the vector 
Vq (see Q35|0 on one hand, and the (largely arbitrary) complex function f(y) characterizing the 
bare solution on the other hand. We note moreover that only in the process with j — 2 there are 
three non vanishing waves in the in-state, a feature which is not found in the one-soliton solution 
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(see below). We also note that not only the deformation F^\y) of the initial profile f(y) (see 
JSSJl) depends on the solution of the ZS equations i|44|) . but also the position of the in-coming and 
out-going solitons (see below) depends on these solutions through their asymptotic behavior, see 
P)a|l . pj)b)l . (|55c|> together with (5HJ, jMbJ, (|54"c|) and HHJ. 

At this point it is evident that numerical experiments, which do not rely on the integration of the 
3WRI system of PDEs Q but just on analytical explicit expressions, are possible if one chooses 
the function f(y) as a "solvable potential", namely a function such that the general solution of the 
ZS equations (|44|) is explicitly known. There exist indeed many such choices. We just pick one 
of them here below to produce explicit formulae which one can use to experiment on three-wave 

interaction proc esses. Before doing this however we observe that the asymptotic profiles Sn (y), 
see ()57all . (|57b(l . (|57cf) and (|59a|) . I|59b|l . (|59c|) . coincide with the asymptotic states of the one- 
soliton solution of the 3WRI equation. This solution obtains by applying the DDT to the trivial 
choice f(y) = and by using the formulae given above. In this particular case the index j which 
we use to distinguish three different initial bare solutions (see is, of course, useless and it is 

dropped. The resulting expression of the one-soliton solution then reads, in our notation, 



\n(x,t) = 2r/g n exp(ipg n x n )x 

x vn~+i 

ljn+l\ 2 ™p(l9n x n ) + \-y n + 2 \ 2 cxp(-?7a„ x n )+\-i n \ 2 cxp [rj (9„ + 2 - 9„+ 1 ) Xn] cxp(-2?7a„ + ig„ + 2 t) ' 



(60) 



where the signs s n , n = 1,2,3, are defined by (|52|l with the order convention (|39|l . The large 
time limits of this solution are readily found to be 

X n(y + c n t,t) — ► 0, (61a) 

t— — S n OQ 



Xn(y + c n t, t) ^ — 2-qg n exp(ipg n y) — ^- — " +2 — f _ , (61b) 



7n+l 7n+2 

t^ s „oo ~^^ ri "' a ' \-f n+1 \ 2 exp(T]g n y) + \j n+2 \ 2 exp(-T]g n y) 

and this last expression coincides, apart from obvious replacement of the parameters 7's, with the 
general expressions (|51a|) with f(y) = or l|51b|) or (|51cl) . If we choose moreover 77 > 0, it is easily 
seen that, because of Q55[l. the asymptotic states here correspond to the expressions given above 
for j = 2 and f(y) = 0. These states may be given also the familiar sech form 

x[~\y) 



_ V9l 


cxp[i (pgi J/+82- 


83)] 




cosh[»7ai(j/-Si)] 




= 0, 






_ V93 


cxp[i (pg 3 y+9i- 


S2)] 




cosh[)j 93(s/-$3)] 




= 0, 








cxp[i (pg 2 y+03- 


0i)] 




cosh [1732(1/— £2)] 





X ( f\y)=0, (62a) 

x[ +) (y) 
xi +) (y) = o, 

by introducing the phases 6 n and the positions £ n via the relation 

r v 1 

7„ = exp(i6>„) exp - (ff n+ i^ n+ i - ff n +2&i+2)J , (63) 
together with the condition 

S16 + 526 + fla& = . (64) 
Moreover, we observe that the asymptotic soliton, whose generic expression is 

7 (+) 7 (-)* 

S(y) = 2r]g exp(ipgy) , (65) 

\T + >\* exp(rjgy) + |-y*. >\ 2 exp{-r)gy) 

may be characterized in the usual way by its amplitude A — E/2, width A = A/E and position 
£ = log(|7'~)/7(+)|) /{rig) where E is the "energy" 

/+00 
dy\S{y)\ 2 = 2\ V g\ . (66) 
-OO 
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Let us now look at one explicit solution. To this aim we choose the undressed wave profile 



f(y) = a exp(-g|y|) , (67) 

where a and q are arbitrary real parameters (a = a* , q = q*) and q is positive, q > 0. The 
chain of steps we have to make to arrive at the asymptotic state expressions 1(5 7a|) . I|57b(l . I(57c(l 
is the following: (i) solving the ZS equations (|44ll and obtaining therefore also the expressions 
of z^^A) and z^\X), see <(46(1 . (ii) computing the vector functions w^\y), see (|53f) and l(45a|) . 
(|45b|) . H45cfl . (iii) computing the constant vector u^^', see (|54a|) . (|54b|) . (|54cl) . and finally (iv) 
using the expressions l|58(l together with (|59a|l . (|59b|l . I(59c() . These steps can be made by using 
the properties of Bessel functions J v (z) (see for instance and we report here only the resulting 
expressions with the following notation: 

C = - , 2 = C cx P(-<7|j/l) , v± = \ ± ^ , (68) 
g 2 2q 

C±=J^ ± (C) - <(C) • (69) 
The solution of the ZS equations 144(1 we choose is 



& fo, A) = exp(-f|y|){i?( 2/ )[C_J^(C)^ + W-C + J_,_(C)^ + W] 

+ ff(-y)[C_J_, + (C)J,_(z) - C+J^COJ-^O*)]} , 

fc(y,A) = cxp(-f|y|){ J ff(y)[C_J, y _(C)^-W + C , + J-,_(C)^-(^)] 

+ i/(-j/)[c_j_, y+ (C)J-„ + (z) + c + J„ + (C)J„ + M]} . 



(70a) 



(70b) 



where i/(y) is the standard Heaviside step function (H(y) = 1 if y > and = if y < 0). 

These expressions of <f)j(y, A) entail that the asymptotic constants (see(@BJ|) may be chosen as 

y (±) i 2 



^=±(-) J„ T (()C T /rK), (71a) 



*2 ±) = (j) J-^(C)G±/r(i/±), (7ib) 

where T denotes the Gamma function (see for instance 6 ). From these expressions all other 
relevant quantities can be computed by simple algebra. For the benefit of the reader we report 
here the standard properties of the Bessel and Gamma functions which have been used to derive 
these expressions: 



' + 0(z 2 ) 



(72) 



zr(^ + )r(^_) [J v+ {z)J v _{z) + J- u+ (z)J_„_(z)] = 2 , (73) 
the first one of which showing the leading term of J v {z) for small \z\. 



It is now rather straight to play with the various parameters, a, q, p, rj, 7„ and c n , together with the 
index j — 1, 2, 3, to explore the interaction of the picked wave profile f(y) = aexp(— q\y\), see (|67|) . 
with solitons in a variety of different regimes with respect to relative velocities, energies and soliton 
positions. Few examples of asymptotic in- and out-state profiles for such processes are graphically 
displayed in the following figures. These examples have been chosen to display processes which 
may be of interest in potential applications of our formulae. They show the effect of collision with 
one soliton (for j = 1, see Fig. 1) and with two solitons (for j — 2, see Fig. 2 and Fig. 3) on 
the exponential picked pulse f(y) l|67|l when such desirable effects as amplitude-amplification and 
width-narrowing of the initial in-state profile f(y) occur. In all figures the modulus is plotted. 

Fig. 1 refers to the case j = 1, see the asymptotic states Q57a|) . with the following values of the 
parameters: c\ = —1, C2 = 0, 03 = 10, p = 0.1, r\ = —1, a = 2, q = 1, 71 = 72 = 73 = 1- In 
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Figure 1: Amplitude amplification. 



particular, Fig. la shows the initial pulse f(y), Fig. lb shows the shape S2 (y), see l(59aj) . of the 
in-state soliton and Fig. lc shows the out-state pulse F^(y) 7 see (158(1 . Fig. 2 and Fig. 3 refer 
instead to the case j = 2, see the asymptotic states (|57bjl . Here two contra-propagating equally 
shaped solitons with speeds c\ = — 1 and C3 = 1 collide with the standing (02 = 0) initial pulse 
f(y). In Fig. 2, with the following values of the parameters: p = 0.1, rj = 1, a = 1, q = 1, 
7i = 72 = 73 = 1, it is shown the case in which the amplitudes and widths of both the two solitons 
and the pulse f(y) coincide, while in Fig. 3, with the following values of the parameters: p = 0.1, 
rj = 2, a = —1, q = 0.6, 71 = 72 = 73 = 1, the amplitude, respectively the width, of the colliding 
solitons is larger, respectively narrower, than that of the initial pulse /(y). 




Figure 2: Width-narrowing and amplitude-amplification. 




Figure 3: Width-narrowing and amplitude-amplification. 



In particular, Fig. 2a and Fig. 3a show the initial pulse f(y), Fig. 2b and Fig. 3b show the 
common shape of the initial solitons s[ 2 \y) and S^\y) (see (I57b|) and (|59b|)) while Fig. 2c and 
Fig. 3c display the out-state profile F^(y), see ()58|l . Figures in the case j = 3 are not reported 
since these collisions feature properties similar to those found in the case j = 1. 

We end this section with the observation that our formulae given here are applicable to several 
other choices of the bare wave f(y) to yield explicit expressions of solutions of the 3WRI equation, 
provided, of course, that the ZS equations ((^|) be exactly solvable. Indeed our choice here, i.e. 
f(y) = aexp(— q\y\), merely indicates how the interested reader may perform simple numerical 
experiments with explicit analytic expressions. 
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4 Conclusion and remarks 



The system is a model equation for the resonant interaction of three waves. It has been 
introduced by the method of multiscale perturbation, as applied to a nonlinear dispersive wave 
equation, with the purpose of describing the time evolution of the amplitudes XiiX2,X3 °f three 
plane waves with wave numbers ki,k2,k% and frequencies u>x, 012,^3 in the resonance conditions 
^1 + ^2 + ^3 = and u>i +u>2 + ^3 = 0. This system is therefore a rather universal model and it finds 
itself in several physical applications. Its integrable version, as characterized by the conditions Q, 
has been therefore actively investigated by means of the theory of solitons. In particular explicit 
solutions have been obtained which correspond to the purely discrete spectrum, i.e. the multi- 
soliton solutions. They may be obtained by repeated application of the DDT (see Section 2) to 
the trivial solution Xi = X2 = X3 = 0. Here we have applied the DDT to a seed solution with 
only two vanishing wave amplitudes while the third one is a non vanishing arbitrary function, 
and we have thereby constructed a larger class of exact solutions of Q outside the class of pure 
discrete spectrum solutions. The feasibility of such construction is due to two special properties 
of the 3WRI system, namely its lack of dispersion and of self-interaction. We have confined our 
construction to the class of localized solutions. However our formulae apply as well to the case 
in which the non vanishing undressed wave f(y), see (|40f) . is outside the class of localized, or 
finite-energy, wave-profiles. The investigation of these "quasi kink" type solutions is left to future 
work. One remarkable example of such solutions has been recently displayed in 7 . The soliton 
solutions introduced there correspond to those obtained by dressing, via the method of Section 
2, the function f(y) = a, a being an arbitrary complex constant ( this case may be viewed also 
as the q — » (singular) limit of our present choice ljfi7)) ). These solutions of the 3WRI equation 
describe the resonant interaction of two "bright" solitons with one "kink" soliton, a process 
which features various behaviours, such as those of boomerons and trappons or the creation and 
annihilation of pairs of bright solitons in the background of a kink. The method of DDT described 
here yields a rather large class of new solutions of the 3WRI equation whose generic behavior is 
of boomeronic type, depending on particular choices of the undressed solution. This is indeed the 
first appearance of boomerons as solutions of a model system of PDEs of such wide applicability 
as the 3WRI equation. In this respect, it should be also pointed out that several systems of 
coupled Nonlinear Schrodinger type equations which possess solutions with similar boomeronic 
phenomenology have been recently found and investigated in [Hj. 
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